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Abstract 

In this paper we investigate qualitative and asymptotic behavior of solutions for a class of 
diffusion-aggregation equations. Most results except the ones in section 3 and 6 concern radial 
solutions. The challenge in the analysis consists of the nonlocal aggregation term as well as the 
degeneracy of the diffusion term which generates compactly supported solutions. The key tools used 
in the paper are maximum-principle type arguments as well as estimates on mass concentration of 
solutions. 



1 Introduction 

In this paper we study solutions of a nonlocal aggregation equation with degenerate diffusion, given by 

Pt = Ap'" + V • (pV(p * V)) in R'' x [0, oo) (1.1) 

with initial data po G L^{W^; (1 + \x\'^)dx) n L°°(R''). Here m > 1, d > 3 and * denotes the convolution 
operator. In the absence of the aggregation term (when y = 0), our equation becomes the well-known 
Porous medium equation (PME): 

Pt - A(p™) = 0. (1.2) 
Note that, formally, the mass of solutions is preserved over time: 

p{-,Q)dx= I p(-,t)dx for aU i > 0. 



Nonlocal aggregation phenomena have been studied in various biological applications such as popu- 
lation dynamics ( |BoCM| . |BuCM| . pM] . [TBL] ) and Patlak-Keller-Segel (PKS) models of chemotaxis 
f |KSj . |LLj . [P]. |FLP| ). In the context of biological aggregation, p represents the population density 
which is locally dispersed by the diffusion term, while V is the interaction kernel that models the long- 
range attraction. Recently, there has been a growing interest in models with degenerate diffusion to 
include over-crowding effects (see for example jTBLj . |BoCMj ). Mathematically, the equation models 
competition between diffusion and nonlocal aggregation. 

In this paper we consider the following two types of potentials: 
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(A) (PKS-modcl) V{x) is a Newtonian potential: 

(1-3) 

1 , , . 

where Cd := 7 ^ — , with ad', the surface area of the sphere § in R . 

(d - 2)ad 

(B) (regularized Newtonian potential) 

V{x)^{J\f*h){x), (1.4) 

where * denotes convolution and h{x) is a radial function in L^(R'* : (1 + \x\'^)dx) fl L°°{M.'^) which is 
continuous and radially decreasing. 

Note that (A)-(B) covers all attractive potentials V with its Laplacian being nonnegative and radially 
decreasing. The restrictions on AV turn out to be necessary for obtaining the preservation of radial 
monotonicity (see Proposition 14. 3p as well as the mass comparison principle in section 5. 

For m > 1, the dynamics in (jl.ip is governed by the energy functional 

Hp) = I i^^r + \p{p * v)\ dx, (1.5) 

and when m = 1, the first term in the integrand is replaced by plogp. Indeed (jl.ip is the gradient 
flow for T with respect to the Wasscrstcin metric (sec for example |AGSj and |CMVj ) . Depending on 
m, the solution of (|l.ip exhibits different behavior. For 1 < m < 2 — 2/d, the problem is supercritical: 
the diffusion is dominant at low concentrations and the aggregation is dominant at high concentration. 
As a result supercritical and critical problems with singular kernels may exhibit finite time blow-up 
phenomena ( [DP] . |HVj . [Sl| . |B1CM] ). On the other hand solutions globally exist with small mass and 
relatively regular initial data, and here the diffusion dominates at large length scale (see [C] and [52]). 
Indeed using the entropy dissipation method f |CJMTU] ) it is shown that the solutions with small 
and norms converge to the self-similar Barenblatt profile f |LSlj - [LS2] and [B2|). 

On the other hand, in the subcritical regime (m > 2 — 2/d), the diffusion is dominant at high 
concentration. For this reason there is a global solution for all mass sizes f|Slj. jBCLj . |BRB| ) . Since 
aggregation dominates in low concentration, one can show that there are compactly supported stationary 
solutions for any mass size (see Proposition 12. ip ). In fact there is no uniqueness result for stationary 
solutions, even for radial solutions, except the well-known result of Lieb and Yau f |LY| ) for the PKS 
model. Furthermore, even for the PKS model, there are few results addressing the qualitative behavior 
of general radial solutions: this is perhaps due to the fact that entropy methods face challenges due 
to the strong aggregation term and the generic presence of the free boundary. This motivates our 
investigation in this paper. 



The main tools in our analysis are various types of comparison principles. While maximum-principle 
type arguments are natural to parabolic PDEs, the classical maximum principle does not hold with 
(jl.ip due to the nonlocal aggregation term, and therefore the standard comparison principle and the 
corresponding viscosity solutions theory do not apply. Instead we establish order-preserving properties 
of several associated quantities: the radial monotonicity (section 4), the mass concentration (section 5), 
and the rearranged mass concentration for non-radial solutions (section 6). 

The following existence and uniqueness results will be used throughout our paper. 
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Theorem 1.1 (Theorem 3 and 7 in jBRBj . Also see |BS| and |S1|). Let V be given by (A) and (B) and 

d > 3. Suppose po is a nonegative function in L^{M.'^; (1 + \x\'^)dx) O L°°{M.'^). Then for m > 2 — 2/d 
there exists a unique, uniformly bounded weak solution p of in x [0,oo) with initial data po. 

Acknowledgments: We thank Scott Armstrong, Jacob Bedrossian and Thomas Laurent for helpful 
discussions and communications. 

1.1 Summary of results 

Let us begin with stating properties of radial, stationary solutions of p.ip : 

Theorem 1.2 (Properties of radial stationary solutions). Let V be given by (A) or (B) and let m > 
2 — |. Let Pa be a non-negative radial stationary solution of (|l.ip with J pA{x)dx ~ A > 0. Then 

(a) Pa is radially decreasing, compactly supported and smooth in its support ( Provosition \2.1\ ): 

(b) PA is uniquely determined for any given A (Theorem 2.2 and Theorem \2.4^ . 

When V is given by (A), the uniqueness of radial stationary solution comes from the well-known 
results of Lieb and Yau ( |LYj ) . Their proof is based on the fact that the mass function satisfies an ODE 
with uniqueness properties; this property fails when V is given by (B). Instead, we look at the dynamic 
equation p.ip . and prove uniqueness out of asymptotic convergence towards a stationary solution. A 
more direct proof of uniqueness and the uniqueness of general (possibly non-radial) stationary solutions 
are interesting open questions. We also mention a recent preprint |BDFj . which studies another type 
of diffusion-aggregation equation: here authors use eigenvalue methods to prove the uniqueness of 
one-dimensional stationary solutions. 

Next we show several results concerning the qualitative behavior of general (nonradial) solutions, 
which will be used in the rest of the paper: 

Theorem 1.3 (Properties of solutions). Suppose m> 1. Let V be given by (A) or (B), and let p{x,t) 
be a weak solution to (jl.ip . which is uniformly bounded in R'^ x [0,T). Then the following holds: 

(a) For any 5 > 0, p is uniformly continuous in R'^ x [S,T); (Theorem \3.1\) 

(b) [Finite propagation property] {p > 0} expands over time period r with maximal rate of Ct^^I"^ 
(Theorem\^; 

(c) If p{-,0) is radial and radially decreasing, then so is p{-,t) for any t £ [0,T) ( Theorem \4. 2^ - 

Both properties (b) (the finite propagation property of the general solutions) and (c) (the preservation 
of radial monotonicity) are new, to the best of the authors' knowledge, for any type of diffusion- 
aggregation equation. For the first-order aggregation equation ( p.ip without the difi^usion term), 
property (c) has been recently shown in [BGL| for the same class of potentials, via the method of 
characteristics. 

We now turn to the discussion of asymptotic behavior of solutions. 
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Theorem 1.4 (Asymptotic behavior: subcritical regime). Let V be given by (A) or (B), m > 2 — |, 
and let p{x,t) be the solution to p.ip with radial, compactly supported initial data Pq{x) G L^{R'^; (1 + 
\x\'^)dx) n L°°(R'^) which has mass A. Let pA be a radial stationary solution with mass A. Then 

(a) The support of p, > 0} stays inside of a large ball {\x\ < R} for allt > 0, where R depends 
on m,d,V and the initial data po f C'orollarv \5.5\) : 

(b) p converges to pA exponentially fast in p-Wasserstein distance for all p > 1 ( Corollarv \5.8\) . and 
\\p{-, t) — pA\\L°°(Rd-j as t ^ oo ( Corollarv \5.9\) . 

The proof of above theorem is based on the mass comparison, i.e. maximum principle arguments 
on the mass concentration of solutions (see Proposition 5.3). The mass comparison property have 
been previously observed for PKS models ( |BKLNj : also see a recent preprint of jCLW| ). However the 
property has not been fully taken advantage of, perhaps because of the success of entropy method for 
the KS model. 

Our method also provides interesting results for asymptotic behavior of radial and non-radial solutions 
in the supercritical regime, when the solution starts from sufficiently less concentrated initial data 
in comparison to a re-scaled stationary profile. (For the definition of "less concentrated than" , see 
Definition 15. ip We point out that in our result the mass does not need to be small as required in 
previous literature (e.g. see |B1| ). and provides an explicit description of solutions which are "sufficiently 
scattered" so that it does not blow up in finite time. 

Theorem 1.5 (Asymptotic behavior: supercritical regime). Let V{x) be given by (A) or (B), and let 
1 < m < 2 — |. Assume po is radially symmetric, compactly supported and has mass A. Then there 
exists a sufficiently small constant S > depending on d, m, A and V , such that if 

Po{X) -< S^pAiSX), 

where " is defined in Definition \5.1\ and paW is given in (|5.35p . then the weak solution p with 
initial data po exists globally and algebraically converges to the Barenblatt profile in rescaled variables 
(Corollary 1^714^ . 

Note that both Thcorcm ll.4l and Theorem ll.5l are limited to radial solutions. The asymptotic behavior 
of non-radial solutions - either in the subcritical regime or in the supercritical regime in terms of the 
re-scaled variables - remains largely open, even for the Newtonian potential. Nevertheless it is possible 
to control the L^-norms of non-radial solutions in terms of radial ones, as we state in the next theorem. 

Let us recall that, for any nonnegative measurable function / that vanishes at infinity, the symmetric 
decreasing rearrangem,ent f* is given by 

/"OO 

r(x) := / Xif>t}'{x)dt, (1.6) 
Jo 

where il* denotes the symmetric rearrangement of a measurable set of finite volume in R'^. 

Theorem 1.6 (Rearrangement comparison and instant rcgularization). Suppose m > 1. Let V be 
given by (A) or (B). Let d > 3 and let p be the weak solution to (jl.ip with initial data poix) G 
L\R'^; (1 -I- \x\^)dx) n L°°(R'^). 
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(a) Let p be the solution to the symmetrized problem, i.e. p is the weak solution to with initial 
data Pq{x). Assume p exists for t G [0,T). Then p*{-,t) -< p{-,t) and \\p{-,t)\\p < \\p{-,t)\\p for 
<t < T and 1 < p < oo (Theorem \6.1\ and Corollarv \6.5\) . 

(b) Suppose m > 2 — 2! then for every < t < 1 we have 

|lp(-,f)Lo.(R.) <c(m,d,A'^^)i"", 
where A = J podx and a d{m-i)+2 ■ ( Pfoposition \6.6\ }. 



Rearrangement results have been obtained before for (|1.2p (Chapter 10 of jV]) and for the two- 
dimensional Keller-Segel model ( |DNRj ) . We largely follow the arguments in [V]. The new component 
in the proof is the introduction of approximate equations to deal with both the degenerate diffusion and 
the nonlocal aggregation term. The L°°-regularization result is interesting on its own: similar results 
have been recently obtained for Kcllcr-Scgel model in }PV) . by a Dc-Giorgi type method. 



2 Properties of the radially symmetric stationary solution 

In this section we consider non- negative radially symmetric stationary solutions of (jl.ip . given by 

^ -p"'-^ + p*V = C in{p>0}, (2.1) 



m — 1 

where we assume m > 2 — |, and the constant C may be different in different positive components of 
p. When V is given by (A) or (B), for any mass ^ > 0, the existence of a stationary solution p with 
mass A is proven in [L] and |B2] . 

Let us define the mass function as follows: 



M{r) := / p{x)dx. 

JB(0,r) 

Since both p and V are radially symmetric, we may slightly abuse the notation and write p * V as a 
function of r. When V = Af, by the divergence theorem and radial symmetry of p and V we have 



f(p*^)M = ^. (2.2) 
or cTdf 



where Cd is the surface area of the sphere S'^ ^ in W^. Similarly, when V is given by (B), for all radially 
symmetric function p, we have that p * ^ is radially symmetric, and 

^ip^V)ir) = ^, (2.3) 
or Odr 

where M(r) := j^^^ r) P* ^^d,x. Note that in both cases, we have dr{p * V^) > 0. 

Proposition 2.1. Let V given by (A) or (B) and suppose m > 2 — |. Then there exists a radially 
symmetric , nonnegative solution p G L^{M.'^) of (|2.ip . Moreover, (a) p is smooth in its positive set; (b) 
p is radially decreasing; and (c) p is compactly supported. 
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Proof. 1. Existence of the stationary solution p follows from [L]: the proof is given in the appendix. 

2. To show (a) for V = M, note that the right hand side of (|2.2p is continuous since /(r) -.^ Jli^'d-i 
is continuous for all r > 0, and /(r) — > as r — > 0. By (|2.2p . p *V is differentiable in the positive set 
of p, which implies that p™^^ (hence p) is also differentiable in the positive set of p. Therefore -^^^ is 
now twice differentiable, hence we can repeat this argument and conclude. When V is given by (B), we 
can apply the same argument on (|2.3p and conclude. 

3. By differentiating (|2.ip we have 

= -^(p*T/) in{p>0}. (2.4) 



m — 1 dr dr 
Due to (|2.2p - (|2.3p the right hand side of (|2.4[) is negative, and thus we conclude (b). 

4. It remains to check (c). Note that (b) yields that p has simply connected support. Hence 
yields 

p(r) = {C - p*V{r))^ . 
When V = Af the proof is similar to that of Theorem 5 in [LYj : since p*V vanishes at infinity, we have 

where Cd is the volume of a ball with radius 1 in M*^. Note that 

p * V{r) < and — p * V{r) ~ ,_ as r — )■ oo. (2-6) 



d-2 
m-1 



If C = 0, implies that 

p(r) = (— p * V{r)) m-i ^ 

where the exponent is greater than — d when to > 2 — |, which contradicts the finite mass property of 
p. Therefore C must be negative and thus p(r) needs to touch zero for some r. 

When V = N'*h,we have p*V = (p^JV) * h. Since h g L^(E'*) and is radially decreasing, using (|2.5p 
we have p * V{x) ~ ^|j_2 as |x| — > oo as well, hence by same argument as above, we can conclude. □ 

Next we state the uniqueness of the radial stationary solution when V = Af. 

Theorem 2.2 ( [LYj ) . Let V — J\f, and suppose m > 2 — |. Then for all choices of mass ^ > 0, 
the radial stationary solution for (jl.ip with mass A is unique. Moreover, the stationary solution is the 
global minimizer for the free energy functional (|1.5p . 



This theorem follows from a minor modification from the proof of Theorem 5 in |LYj , which proves 
uniqueness of the stationary solution of a slightly different problem. Their proof consists of two steps: 
firstly for a given mass, they first show the global minimizer of (jl.Sp is unique, and secondly they prove 
every radial stationary solution is a global minimizer for some mass. Theorem l2.2l and the homogeneity 
of N yields the following: 
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m > 2 



m = 2 



2-§ < m < 2 



Figure 1: Stationary solutions with different mass for different to, wliere J pAdx < J psdx. 

Corollary 2.3. Let V and to be as in Theorem \2.2\ and let pM be the radial solution of (|2.1|) with 
mass M . Then 

Pm{x) = api{a ~x) with a := M <i{rn-2+2/d)) ^ (2.7) 

In particular if A < B then maxpg > maxp^ and the following dichotomy of behavior is observed.: 
(see Figure 1 ). 

(a) When m > 2, {pA > 0} C {p^ > 0}. 

(b) When 2 - I < TO < 2, {pe > 0} C {p^ > 0}. 

Wc point out that the proof in jLY| cannot be generalized when V is given by (B): the difficulty 
lies in the second step. When V = M, for any radial stationary solution p, its mass function M{r) := 
I\x\<r p{x)dx solves a second order ODE 

/ TO , M'{r) yn-i\' _ M(r) 

where M(0) = is prescribed. It follows that A/(r) is unique for a given p{Q) ~ linv-^o M' [r] / {adr"^^^), 
which implies that p can be uniquely determined by p{Q). This property then allows both the radial 
stationary solutions and the global minimizers to be parametrized by their values at the center of mass 
(see Lemma 12, |LY| ) . 

When V is given by (B), M{r) solves a nonlocal ODE, hence different stationary solutions may have 
the same center density: thus the above argument in |LYj cannot be applied to prove the second step, 
necessitating an alternative approach. Instead of dealing with the stationary equation (j2.ip directly, we 
will consider the dynamic equation (jl.l[) and prove the uniqueness of the radial stationary solution by 
their asymptotic convergence. Indeed the following theorem is one of the main results in our paper. 

Theorem 2.4. [Corollary \5.10f Let V be given by (B), and suppose to > 2 — |. Then for any A> Q, 
the radial stationary solution of (jl.ip with mass A is unique. 



3 Qualitative properties of solutions 

In this section several regularity properties, including the finite propagation property, will be derived for 
general weak solutions of (jl.ip . We point out that the results in this section hold for general (non-radial) 
solutions. 
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Theorem 3.1. Suppose m> 1. Let V be given by (A) or (B), and let p be a weak solution of (1-1) with 
its initial data pQ bounded with compact support. Further suppose p is uniformly bounded in M.'^ x [0,r]. 
Then 

(a) For any S > 0, p is uniformly continuous in x {d,T]. 

(b) [Finite propagation property] Suppose {x : p{-,t) > 0} C Br{0). Then 

{x : p{-, t + h)>0}c Bji+chU2 (0) forO<h<l, 
where the constant C > depends on m, d, po and ||AV^||i. 

Proof. 1. Let us consider the case V = TV. This is the most singular case and paraUel (and easier) 
arguments hold for V given by (B). Let 

Co = snp{p{x,t) : {x,t) G M" x [0,T)}. 

Observe that treating the convolution term $ := y * p as a priori given, p solves 

A(p™) + V-(pV$). (3.1) 

Also, for aU t e [0,r), $ satisfies 

|V$|(-,i) <Co / \V^^\iy)dy+\\p{■,t)\\l sup |VA/-(y)| < Ci, (3.2) 
"'|y|<i |y|>i 

where Ci depends on Co, the and sup- norm of p, and the dimension d. Also 

|A$|(-,t) < IIpIIloo < Co for aU t e [0,r). (3.3) 

The bounds p.2|) - p.3|) and Theorem 6.1 of |Dibj yields the uniform continuity of p in M'' x [5,T). 

2. Next we prove (b). First of all let us point out that the standard comparison principle holds between 
weak sub- and supersolutions of (|3.ip . For the case of time-independent potential ^{x,t) = $(a;). 
Proposition 3.4 of jBHj asserts that the comparison principle between weak sub- and supersolution of 
p.l[) holds if the potential function $ is independent of the time variable and | V$|, | A$| < C. The proof 
in |BH] is based on an approximation of the original problem p.ip by a sequence of regularized problems 
which satisfy the comparison principle (see sections 4, 7, 8 of [BHj ). This argument straightforwardly 
extends to our (time-dependent potential) case, and one can verify that comparison principle between 
weak sub- and supersolution of (|3.ip holds. 

We will now construct a supersolution of p.ip to compare with p over a small time period to prove 
the finite propagation property. First observe that the pressure function defined by u := ;;friP"'^^ 
formally satisfies the PDE 

ut = {m- l)uAu + |Vup + Vu • + (m - 1)mV$. (3.4) 

Based on this observation, we will first construct a supersolution of (|3.4p , and use the pressure-density 
transformation to construct the corresponding supersolution of (j3.ip . Let us define 

U{x,t):^A inf e'^Wyl + cut - B)., 

\x-y\<C-Ct 
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where w = (1 + (m — l)(d — I)) A, and the constants B and C wiU be chosen later. 

Let E :~ {\x\ < 2B} x [0,uj^^B]. Due to Proposition 2.13 in jKL| . [/ is a viscosity (or weak) 
supersolution of (|3.4p if C is chosen to be larger than max(Co, Ci) given in p.2p - (|3.3p . In other words, 
U satisfies 

t7( > (m - 1)UAU + I V?7p + C\\7U\ + CU in {[/ > 0} n S, 
and the outward normal velocity Vx.t of the set {U > 0} at {x,t) G d{U > 0} satisfies 

V^^t^to + C >A + C > |Vf/| + C. 
Hence p := (iri^f/)i/(™-i) satisfies 

/5t > A(p'") + qvpi + -^p 

771—1 

in the domain E, in the viscosity sense (see jKLj for the definition of viscosity solutions of ((SU])). 

Moreover, observe that U is Lipschitz continuous in space, and continuous in space and time. 

Using this regularity of p as well as the above estimates on the derivatives of $, it follows that p is a 
weak supersolution of (|l.ip in S, if we choose C greater than (m — l)Ci. More precisely the following 
is true: for all times < < < lo^^B and for any smooth, nonnegative function 
ipix.t) : K" X (0,cx}) ^ M with {i/'(-,t) > 0} C {\x\ < 2B] for < < < uj'^B, we have 

p(-, t)^{-, t)dx > J p{-, 0)7/'(-, 0)dx + j j iP'^'M^ + Pi't - /5V$ • \/4')dxdt. 



Now suppose {p{-,to) > 0} C Br{0) for some to g [0,r]. Let us compare p with p in 
E := {|a;| < 2B} x [to,to + h], with B = R + h^^^ and A = 2Coh-^''^. 

Since {p(-,to) > 0} C Bii{0) and p < Co, we have p < p at t — to- Therefore comparison principle 
asserts that p < p a.t t = to + h. In particular 

{pi-, to + h)>0}(l {p{; to + h)>0} = Bj,^^,y^./. (0), 

where M = h^/'^ + /iw and a; = (1 + (m - l){d - l))yl. This proves (b). □ 

Remark 3.2. Due to \BRB^ . when 771 > 2 — ^, there exists a global weak solution p of (1.1) with initial 
data po- Moreover, p is uniformly bounded in X (0,00) due to Theorem 10 in JBRBf . so in that case 
we may let T = 00. 



We finish this section with an approximation lemma which links case (A) and (B). Let 

■.= e-^hi-) 
e 

with h being the standard moUifier in M'' with unit mass, and let be the corresponding solution 
of (|1.1|) with V = JV * h"^ and with initial data po- Then Lemma 8 in |BRBj yields that {p'}e>o are 
uniformly bounded for t £ [0, T] for some T. This bound as well as Theorem 6.1 of |Dibj yields that the 
family of solutions {p'^} are equi-continuous in space and time. This immediately yields the following 
result: 

Proposition 3.3. Let po be as given in Theorem \3.1\ Let V = Af * h'^ and let p^ be the corresponding 
weak solution of (jl.ip with initial data po- Let p be the unique solution to (jl.ip with V = M and initial 
data Po, and assume p exists for t £ [0,T), where T > may be infinite. Then the solutions p"^ locally 
uniformly converge to p mR'^ X [0,T). 
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4 Monotonicity-preserving properties of solutions 



In this section, we show that when V is given by (A) or (B), solutions with radially decreasing initial 
data remains radially decreasing for all future times. The main step in the proof is the maximum 
principle-type argument applied to the double-variable function 

:= p{x,t) - p{y,t) in {|a;| > \y\} x [0,oo) 

to ensure that ^ cannot achieve a positive maximum at a positive time. 

We begin with an observation on the convolution term; the proof is in the appendix. 

Lemma 4.1. Let V{x) be given by (B). Let u{x) be a bounded non-negative radially symmetric function 
in with compact support. Further suppose u{x) is not radially decreasing, i.e. there exists Oi = 
(a, 0, 0) and bi = {/3, 0, 0) with a, f3 > such that 

— u(ai) = sup u(6) — u(a) > 0. (4-1) 

\a\<\b\ 

Then we have 

{u^AV){b,)-{u*AV)iai) < \\AV\\L^iu{b,)~u{a,)). 
Theorem 4.2. Let V{x) be given by (A) or (B). Suppose that the initial data 

p{x, 0) S L^{M.'^; (1 -I- |a;p)(i2:) fl L°°(R'^) is radially decreasing, i.e. p{x, 0) is radially symmetric and is 
a decreasing function of\x\. We assume a weak solution p exists for t E [0, T), where T may be infinite. 
Then p{x,t) is radially decreasing for all t G [0,T). 

Proof. 1. Without loss of generality we assume that V is given by (B), and p{x,0) is positive and 
smooth. Then a classical solution p{-,t) exists for all t > 0, and we want to show p{-,t) is radially 
decreasing for alH > 0. When V = JV, we can use mollified Newtonian kernel to approximate A/"; and 
for general radially decreasing initial data, we can use positive and smooth functions to approximate 
p{x,0). Then the result follows via Proposition 13.31 

2. Radial symmetry of p for all i > directly follows from the uniqueness of weak solution. To prove 
that p is radially decreasing for all time, let us define 

w{t) := sup p{b,t) — p{a,t). 

\a\<\b\ 

Since p is uniformly bounded and uniformly continuous in x [0,oo), w{t) is continuous in t, and 
uniformly bounded for t S [0, oo). Moreover, p(x, 0) being radially decreasing is equivalent with w{Q) ~ 
0. We will use a maximum principle-type argument to show that w{t) = for all t > 0, which proves 
the theorem. 

Suppose w ^ 0. Then for any A > the function w(i)e~^* has a positive maximum at t = ti for 
some ti > Q. We will show that this cannot happen when we choose A > 2||p||oo|| Ay||i. 

At t ~ti there exists ai = (a, 0, 0) and bi = {13, 0, 0) such that a < /3 and 

p{b,,t,)- p{ai,t,)=w{ti)>0. (See Figured) (4.2) 
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ai bi 
Figure 2: Graph of p at time ti 



Moreover by definition p{bi,t) — p{ai,t) < w{t), and thus 

-/9(ai,t))e-^*) =0 at t = 

which means 

ptibi,ti) - pt{ai,ti) = X{p{bi,ti) - p{ai,ti)). 

Further observe that has a local minimum (in space only) at ai and a local maximum 

This yields 

Vp(ai,ti) = a7T,dVp{bi,ti) = 0, 

as well as 

Ap'"(ai,ii) > and Ap'"(6i,<i) < 0. 

Let us now make use of the equation (1.1) that p satisfies to get a contradiction: we have 

Ptibi,h) - Ptiai,h) = Ap"(6i,ti) + V-(pV(p*F))(6i,<i) 

-Ap™(ai,ti) - V • {pW{p*V)){ai,h) 

< p{biM){p * ^V){biM) ~ p{aiM){fi * Al/)(ai, ti) 

= p{biM)[{p * ^V){bi,h) - (P * ^V){auti)] 
+ {p{bi,ti) ~ p{ai,h)){p^ AV){ai,h). 

In order to bound the first term in (j4.4p . we apply Lemma l4.11 which gives 

{p * AV){bi,h) - {p * Ay)(ai, ii) < II Ay ||i(p(6i, ii) - p(ai, ti)), 
and for the second term we use 

(p*AV){ai,h) < ||p|U||Ay||i. 
Due to the estimates dm yields that 

Pt{bi,h)-pt{ai,h) < 2\\p\\^\\AV\\i{p{bi,ti)- p{ai,h)), 
which contradicts (|4.3p since we chose A to be strictly greater than 2||p||oo|| Ay||i. 
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The following proposition states that in the previous theorem, the eondition that AV is radially 
decreasing is indeed neeessary. 

Proposition 4.3. Let V{x) be radially symmetric, AV > 0, with AV continuous, but not radially 
decreasing. Then there exists a radially decreasing initial data Pq such that the solution p{x,t) of (1-1) 
starting with initial data po does not preserve the radial monotonicity over time. 

Proof. Since AT^ is not radially decreasing, we can find xi,X2 G K'^, such that < \xi\ < \x2\; and 
AV{xi) < AV{X2). 

For a small e > 0, let po{x) be given as below: 

1 X 

Po{x) = eXB{o,x2+i) * (l>{x) + -^M-)-^ 

where xe is the characteristic function of E and is a radially symmetric mollifier with unit mass 
and supported in 5(0, ro), where rg < min{l, |xi|/2}. Note that in a small space-time neighborhood of 
xi and X2, p solves a uniformly parabolic equation, and thus is smooth. 

Since Ap"'{x^, 0) = \7p{x^,0) = for i = 1, 2, we have 

pt(xi,0) = po{xi){po * AV){x^), i = 1,2. 

Since po{xi) = po{x2), to show pt{xi,0) < pt{x2,0), it suffices to prove 

{po*AV)ixi)<{po*AV){x2). (4.7) 

Note that pQ * AV locally uniformly converges to AV{x) as e ^ 0. Since AV{xi) < AV{x2), if we let e 
be sufficiently small, we would have (|4.7|) . In particular p{xi,t) < p{x2,t) for small t > 0, which means 
p{x, t) stops being radially monotone as soon as t > 0. □ 

5 Mass Comparison and asymptotic behavior for radial solu- 
tions 

Recall that there is no classical comparison principle for (1.1), due to the nonlocal term; however, we 
will prove that a comparison principle actually hold for the mass function 

M{r,t) = M{r,t;p) := / p{x,t)dt (5.1) 

Js(0,r) 

(see Proposition [531) • We point out that the corresponding property has been observed for p.2p (jV]) 
and also for the Keller-Segel model f |BKLNj ). It turns out that mass comparison holds for (|l.ip if 
the potential V satisfies AV > in the distribution sense (Proposition 15.31 ) As we will see later, 
mass comparison effectively describes the asymptotic behavior of radial solutions in both sub- and 
supercritical regime. 
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5.1 Mass comparison 



First note that, if p is a weak solution of in the positive set, p is a bounded solution of a locally 

uniformly parabolic, divergence- type equation with continuous coefficients. Hence due to [LSUj . p is at 
least in space and time variables. It follows that M{r,t) is in space and in time in {p > 0}. 

Let us denote the support of p at time t by B{0, R{t)). Then we compute the PDE that M{r,t) 
satisfies in {?- < R{t)}. Due to p.3p . we have 



dt 



'dB{0,r} 

^_^rd, dM 1 , dM 1 w ^ 



where 

M{r,t)=M{r,t;p):= I {p{- ,t) * AV){x)dx. (5.3) 

"'B(0,r) 

Definition 5.1. Let pi and p2 be two non-negative radially symmetric functions in L^(R''). We say 
pi is less concentrated than p2, or pi -< p2, if 



pi{x)dx < / p2{x)dx for all r > 0. 

B(0,r) JB(0,r) 

Definition 5.2. Let pi{x,t) be a non-negative, radially symmetric function in L°°{{[0,T]; L^{M.'^) D 
L°°{M.'^)), which is in its positive set. We say pi is a supcrsolution of (j5.2[) if Mi{r,t) :~ M{r,t; pi) 
and Ml := M(r,t; pi) satisfy 

dM, ^_^d/,dMi 1 .,n\ ,dMi 1 

in the positive set of pi . 

Similarly we define a subsolution of (|5.2p . 



Proposition 5.3 (Mass comparison). Suppose m > 1. Let V be given by (A) or (B), and let pi{x,t) 
be a supcrsolution and p2{x,t) be a subsolution of ()5.2p for t G [0,r]. Further assume that the pi's 
preserve their mass over time, i.e., J pi{-,t)dx and J p2{-,t)dx stay constant for all <t <T. Then 
the mass functions are ordered for all time: i.e., if p2{x,0) -< pi{x,0), then we have 

P2{x, t) -< pi{x, t) for all t G [0, T]. 



Proof. Wc claim that M2{r,t) < Mi{r,t) for all r > 0,t G [0,T], which proves the proposition. At 
t = 0, we have M2{r, 0) < Mi{r, 0) for all r > 0. For the boundary conditions of Mi, note that 

Mi(0, t) = A/2(0, t)=0 for aU t G [0, T]; 

lim,^^(Mi(r,0 ~ M2{r,t)) = J^Mx,0) - p2{x,Q))dx > for aU t G [0,T]. 
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To prove the claim, for given A > 0, we define 

w{r,t) := {M2{r,t)-Mi{r,t))e-^\ 

where A is a large constant to be determined later. If the claim is false, then w attains a positive 
maximum at some point (ri, ti) £ (0, oo) x (0, T]. Moreover, since the total masses of pi's are preserved 
over time and thus are ordered, we know that (ri,ti) must lie inside the positive set for both pi and 
P2, where the Mi's are Cl']. 

Since it is assumed that w attains a maximum at (ri, ti), the following inequalities hold at this point: 

d(M2 " Ml) 

wt>0 =^ ^—Ql > - Ml); (5.5) 

dMi dM2 , , 

-^ = 0=^^-^; (5.6) 

d'^Mi a^Ms , , 

Now consider the first term on the right hand side of (|5.4p . and the corresponding inequality for M2. 
jS^ and dSTT]) imply that 

d (.dMi 1 „\ d f.dM2 1 ,,„\ w \ o\ 
Subtracting the inequality (j5.4p with the corresponding inequality for A/2, and using (|5.8p . we obtain 



We next claim 

{M2 - Mi)iri,ti) < C{M2-Mi){ri,ti), 

where C only depends on V. For the Newtonian potential this is obvious, since M = M. For a mollified 
Newtonian potential, we estimate (A/2 — A/i)(ri,ti): 



A/2(ri,ti)-A/i(ri,ii) = / {{p2 - pi) * AV) XB(o,r,)dx 

{P2 - Pl)(XS(0,ri) * l^V)dx. 



By our assumption, W is radially decreasing, thus x_B(o,ri) is radially decreasing and has maximum 
less than or equal to || Ay||i. Therefore we can use a sum of bump functions to approximate the function 
XB(o.ri) * Ay 1 where the sum of the height is less than || Ay||i. Hence, like in the proof of Lemma WA[ 
we get 

M2{riM) - Mi{ri,ti) < I|Ay||i sup(A/2 - A/i)(x, ti) = I|Ay||i(A/2 - A/i)(ri, ii). (5.10) 

X 

We plug the estimate (|5.10p into (|5.9|) . which then becomes 

^^^^i^^— ^<Pi||AT/I|i(M2-Mi) at(ri,ii). (5.11) 
at 

Now choose A > supjj„^(o,T] II^^IIiII/Oi|Il°°(e<'x[o,t])- Then we see that (|5.1ip contradicts (|5.5p . □ 
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Observe that p.ip can be written as a transport equation 

pt + V ■ {pv) ^ 0, 

where the velocity field v is defined by 



vix, t- p) := -V(p™-i) - V(p * V). (5.12) 

m — 1 



Then the mass function for a radial solution of (jl.ip satisfies 



M(r,0 = -p(r,i) / (« -71)^5. (5.13) 

JdB{0.r) 

The above observation along with Proposition 15.31 immediately yields the following corollary: 

Corollary 5.4. Suppose m > 1. Let V be given by (A) or (B). Let po{x) be a continuous radially 
symmetric function, which is differentiable in its positive set. We assume that the velocity field of po is 
pointing inside everywhere, i.e., for v as defined in (|5.12p . 

z/(x;po)- (-^) >0 m{po>0}. (5.14) 
Let p be the weak solution of (|l.ip with initial data po -< p(-,0). Then po -< p{-,t) for all t>Q. 
Proof. Let us define 

pi{x,t) := po{x) for {x,t) £ R'^ x [0,oo). 

Then (|5.13p and (|5.14[) yield that pi is a subsolution of (|5.2p . Therefore, Proposition 15.31 applies to p 
and pi and so we are done. □ 

As an application of Corollarv l5.4[ we will show that when the initial data is radially symmetric and 
compactly supported, the support of the solution will stay in a large ball for all times. 

Corollary 5.5 (Compact solutions stay compact). Let V be given by (A) or (B), and /ei m > 2 — |. 
Let p solve (|l.ip with a continuous, radially symmetric and compactly supported initial data p(a;, 0). 
Then there exists R> depending on m,d, ||Ay||i and p(-,0), such that 

{p{-,t) > 0} C {\x\ < R} for all t > 0. 

Proof. 1. We will first assume that p(0, 0) > 0. Let A J-g^^ p{x,Q)dx, and let pa{x) be a radial 
stationary solution with mass A. For any continuous radial initial data with p(0, 0) > 0, we can choose 
a > sufficiently small, such that 

pi{x,t) := a'^pAiax) -< p{x,0). 
Our aim is to show that the velocity field of pi{x,t) is pointing towards the inside all the time, i.e., 

v{r,t;p{) := v{r,t-p,) ■ ^ = |:Pr'M + ^^{pi*V)> 0. (5.15) 
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Let us assume that V is given by (B); the argument for V given by (A) is parallel and easier. Recall 
that the stationary solution pAix, t) satisfies the following equation in its positive set: 



> 



TO — 1 dr (TdT'^ ^ 

Therefore it follows that 

™ ^ n^-^(r.\ - Ain-i)d+i (^Pa ^ _ „(m-i)d+i ^^(°-^;P^) 

m — 1 or m — 1 ar (Td(arj" ^ 

Secondly observe that M{r,t;pi) satisfies 

M{r,t;pi) = / a'^pA{ay)AV{y-x)dydx 

J B{0,r) JM'' 

I pA{y)a-'^AV{a-\y-x))dydx 

{0,ar) JR<1 

/ PA* AVdx (since a^'*Ay(a~^x) >- when < a < 1) 

JB{0,ar) 

= M{ar;pA). 

(Note that when V is given by (A), direct computation yields AI{r,t;pi) = M{ar; pa)-) 
Due to (|2.3p and the above inequalities, it follows that 

<r,t;p,) = |-p™-i(,) + |-(p,*V) 

> prHr) + ^^ (5.17) 

> (l-a^(™-2+2/d))^d-iM(ar^_ 

Since to > 2 — 2/(i, the above inequality yields that the inward velocity field v{r,pi) > when a < 1. 
Therefore Corollary 15 .41 implies that p{-,t) >- pi for all t >0. Since p and pi have the same mass A, it 
follows that 

{p{-,t) > 0} C {pi(-,0 > 0} for all t > 0, 

and we can conclude. 

2. The assumption p{0,0) > can indeed be removed, since p{0,t) would still become positive in 
finite time even if p{0,0) ~ 0. This is because, for (|1.2p . it is a well-known fact that the solution will 
have a positive center density after finite time: this can be verified, for example, by maximum-principle 
type arguments using translations of Barenblatt solutions. 

Note that a solution of (|1.2[) is a subsolution in the mass comparison sense. Hence one can compare 
p with a solution of (|1.2p with initial data po and apply Proposition 15 . 31 to conclude that ip ~< p. Now 
our assertion follows due to the continuity of and p at the origin. □ 
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5.2 Exponential convergence towards the stationary solution in the subcrit- 
ical regime 



As an application of Proposition 15.31 we will prove the asymptotic convergence of radial solutions to 
the unique radial stationary solution when the potential is given by (A) or (B). 

Theorem 5.6 (Exponential convergence of radial solutions with the Newtonian potential). Let m > 
2 — I and let V be given by (A) or (B). For given po > 0; a continuous, radially symmetric function 
with compact support, let p{x,t) be the solution to (|l.ip with initial data pQ. Next let pAix) be a radial 
stationary solution with mass A := J po(x)dx. Then M(r,t) := M(r,t;p) satisfies 

\M{r,t)~M{r-pA)\<Cie-^\ 

where Ci depends on po, A, m, d, V , and the rate A only depends on A, m, d, V . 

Proof. 1. We will only prove the case when V satisfies (B); the case for (A) can be proven with a 
parallel (and easier) argument. Note that we may assume po{0) > since otherwise p{0,t) will become 
positive in finite time as explained in step 3 of the proof of Corollary 15.51 

2. Let Pa be a stationary solution with the same mass as po, given as in the proof of Corollary 
15.51 Since po is compactly supported, continuous and with po(0) > 0, we can find a sufficiently small 
constant a > such that 

a'^PA{ax) -< Po and a^'*p^(a^^a;) >- po. 

3. With the above choice of a, we next construct a self-similar subsolution (j){x,t) of (|5.2p with initial 
data 4>{x, 0) = a'^pA{ax) such that M^{-, t) := M{-, t; (f) converges exponentially to A/(-; pa) as i -> oo. 

Here is the strategy of construction of 4>{x,t) . Due to (|5.18|) . for all < a < 1, the inward velocity 
field v{r) := v{r;a'''pA{ax)) given by (|5.15p satisfies 

v{r) > (1 - a^(m-2+2/d)y^ M{ar;pA) ^ ^ 

dM{ar;pA) 

Observe that — — equals the average of pA * Ay in the ball { x < ar\. By Proposition 12. 1[ pA 

ad[arY 

(hence pA * Ay) is radially decreasing, and thus we have 

C,<^^^^p^<C, in{p^>0}, (5.19) 

where Ci, C2 only depend on A, d, to, V. This gives a lower bound for the inward velocity field v 

v{r) > Cia'^il - a''(™-2+2/<i)^^._ ^5 20) 

We will use the above estimate to construct a subsolution 4){r,t) of (|5.2p . Let us define 

<i>{r,t)^k''{i) pA{k{i)r), (5.21) 
where the scaling factor kit) solves the following ODE with initial data fc(0) = a: 

k'{t) = Ci(fc(i))''+i(l - (fc(t))''(™-2+2/d)-)^ ^5 22) 
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Since m > 2 — 2/ d, k'{t) > when < fc < 1, and since fc = 1 is the only non-zero stationary point for 
the ODE for < fc(0) < 1 we have Umt^oo k{t) = 1. Since 

Clk'^il - /t'*("-2+2/d)) ^ _Cld(771 - 2 + 2/d){l -k) + o(l - fc), 

it follows that 

< 1 - k{t) < e-Cid(™-2+2/d)t^ (5 23) 

which implies 

< M(r; pa) - M^{r, t) < e-Ci'i('"-2+2/d)t_ (5 24) 

Next we claim that cj) is a. subsolution of (|5.2p . i.e., 

1 d /',dMj, 1 fdMj, 1 \ ~ , , , , 



To prove the claim, first note that by definition of (j){r, t) we have M^{r, t) = M {k{t)r; pa)- Hence, due 
to (|5.22p and definition of (j>, the left hand side of (|5.25p can be written as 



r 

^ir,t) = drMik{t)r;pA) k'{t)r 

■•^PA{k{t)r)k-^^^{t)k'{t) = adrU{r,t)Cik'^{l - k'^i"^-^+Vd)y (5 25) 



DM, 
~dt 



On the other hand, we can proceed in the same way as (|5.20p . replacing a by fc, to obtain 



Therefore 



m — 1 dr ad.r'^ 



RHSof(I03 = cFcir'^^^-^^ 

or 



m — 1 dr (JdT'^ ^ 

> arf/0Cifc''(l - A;''("-2+2/d))^ 

thus indeed satisfies (|5.25p . and the claim is proved. 

4. Similarly one can construct a supersolution of (|5.2p . Let us define 

r7(r,t) ■.= k\t) pA{k{t)r), 
where k{t) solves the following ODE with initial data fc(0) = i; 

k'{t) = C2fc-''+^(l - fc'i('"-2+2/d))^ 

where C2 is defined in (|5.19p . Arguing parallel to those in as in step 3 yields that 77 is a supersolution 
of (|5?2l) and 

< Mr,{r, t) ~ M{r; pa) < e-C2'i('"-2+2/<i)t^ foj. all r > 0. (5.27) 

5. Lastly we compare ^,77 with the weak solution p of Since 

</)(•, 0) -< p{-,0) -< T]{-,0) (see Figure E]), 
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Figure 3: Initial data for cf), p and rj 



Proposition 15.31 yields that 

M^{-,t) < M{;t) < M„(-,t). (5.28) 

By ([5?24| and ([5?27| . wo obtain 

\M{r,t) ~ M{r;pA)\ < e-<^i''('"-2+2/d)t 

□ 

Using the explicit subsolution and supersolution constructed in the proof of Theorem 15.61 we get 
exponential convergence of p/A towards pa/A in the p-Wasserstein metric, which is defined below. 
Note that the Wasserstein metric is natural for this problem, since as pointed out in [AGSj and jCMVj , 
the equation is a gradient flow of the energy (|1.5|) with respect to the 2- Wasserstein metric. 

Definition 5.7. Let pi and p2 be two (Borel) probability measure on with finite p-th moment. Then 
the p- Wasserstein distance between pi and p2 is defined as 



Wp{pi,p2) --^ ( inf 1/ \x - y\PTT{dxdy)\) 



where V{pi,p2) is the set of all probability measures on R'^ x K.'' with first marginal pi and second 
marginal p2 ■ 

Corollary 5.8. Let p, pA, A,Ci, X be as given in Theorem \5.b\ Then for all p > 1, we have 



At 



Proof. Before proving the corollary, we state some properties for Wasserstein distance, which can be 
found in jv!]. For two probability densities /o, /i on R'', the p- Wasserstein distance between them 
coincides with the solution of Monge's optimal mass transportation problem. Namely, 

r — 

Wp{fiJo)^(ini fo{x)\x-T{xWdxy, (5.29) 
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where T is a map from to M'', and T^fo = fi stands for "the map T transports /o onto /i", in the 
sense that for ah bounded continuous function h on R'', 

/ h{x)Ji{x)dx= I h{T{x))fo{x)dx. 

Let (f> be the subsolution constructed in the proof of Theorem 15. 6[ where we proved the radius of 
supp (j) converges to the radius of supp pA exponentiaUy in time. Note that </)(-, t) is a rescahng of pA, 
hence the convergence of support imphes that there is a map T^[-,t) transporting (f>{-,t) onto pA with 
sup3,g{p(. f-j^Qj. \x — Ttj,(x)\ decaying exponentiaUy in time. Once we find such T^, we can use (|5.29p to 
show that Wp{(j){-,t), pa) decays exponentially. 

Without loss of generality, we assume the mass A = 1 for the rest of the proof to avoid dividing p 
and Pa by A every time. The transport map can be explicitly constructed as 

T4x,t) = -^^M^\M{\x\;pA),t). 

Recall that (/) is a rescaling of pA, defined as (j){x,t) = {k{t)Y pA[k[t)x), where k{t) < 1 for all t, 
and k{t) converges exponentially to 1. In this case the defined above can be greatly simplified as 
T^{x,t) = x/k{t). That gives us the following upper bound bound for Wp{4>{- ,t), pa)'- 

where R is the radius of support of pA- Due to the estimate of k{t) in (j5.23p . we obtain the exponential 
decay 

Wp{(t>{-,t),PA)<Cie-^K 
We can apply the same argument to the supersolution 'q[-, t) as well. 

To show that Wp{p{- ,t), pa) decays with the same rate, it is natural to consider the following map 
T{-,t) which transports pA onto p{-,t): 

T{x,t) = ^M-\M{\x\-pA),t). 
\x\ 

Then we have \T{x,t) - x\ = \M-^{M{r; pA),t) - r\, where r = |a:|. Due to ((5?28l) . we have 
which gives 

\T(x, t) ^ x\ < max{|T0(a;, t) ~ x\, |T,,(x, t) — a;|}. 

Hence we conclude that 

WpipC^t)^PA) < Wpi4>{;t),pA)+Wp{r^{;t),pA) < Cie"^*. 

□ 

In fact one can also obtain the uniform convergence of p(-, t) to pA in sup-norm, however the conver- 
gence rate would depend on the modulus of continuity of p. Theorem 15.61 and the uniform continuity 
of p and Pa, as well as the fact that pA is compactly supported, yield the following: 
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Corollary 5.9. Let p,pA, Ci and A he as given in Theorem \5.6[ Then we have 

lim ||p(a;,i) - pA(a;)||i,~(Rd) = 0. 

Note that uniqueness of pA is not required in the proof of Theorem 15.61 Indeed, uniqueness of pA 
follows from the asymptotic convergence of p: if there are two radial stationary solutions p\ and p\ 
with the same mass, Corollary 15.91 implies p{-,t) — ^ in L°° norm for i = 1,2 when p is given as in 
Theorem 15.61 This immediately establishes the uniqueness of the radial stationary solution. 

Corollary 5.10. Let V be given by (A) or (B), and /et m > 2 — |. Then for all A > 0, the radial 
stationary solution pA for with J pA{x)dx = A is unique. 



5.3 Algebraic convergence towards Barenblatt profile in the supercritical 
regime 

In this section, wc consider the asymptotic behavior of radial solutions in the supercritical regime, i.e. 
for 1 < 771 < 2 — I . In this case the diffusion overrides the aggregation and thus the solution is expected 
to behave similar to that of Porous Medium Equation (PME) in the long run. In fact recently it was 
shown in |B1| (and also in |Slj ). by making use of entropy method as well as functional inequalities, 
that the solution of (1.1) with a general class of V and with small mass and small iji'^-™-)^/'^ norm 
converges to the self-similar Barenblatt solution U{x,t) with algebraic rate, 

U{x,t) = t-f^'^iC - (5.30) 

2777, 

where C is some constant such that ||Z//(-,0)||i = ||/r7(-, 0) 



Here we will give a complementary result to jBl| and |Slj in the case of radial solutions, by using 
mass comparison (Proposition 15. 3p . We point out that in our result the mass does not need to be small, 
and we provide an explicit description of solutions which are "sufficiently scattered" so that they do 
not blow up in finite time. Of course the method presented in |B1| is much more delicate and yields 
optimal convergence results for general solutions with small mass in the supercritical regime. 

Let p be the weak solution to (|l.ip . Following jV], we re-scale p as follows: 

p{X,t) = {t+l)''p{x,t); X = x{t + 1)-I^: r = ln(t + l), (5.31) 

where 

P = a/d. 



d{m - 1) + 2 

Then /i(A, 0) = p(x, 0), and /i(A, r) is a weak solution of 

lir = A/i™ + /3V ■ (mV^) + e(i-")"V • {^lV{p * (A/- * h{\ r))), (5.32) 

where 

h{\T)-e'^P^l^V{\e^^). (5.33) 
(When V ~ N one should replace the last term by e'^^~"^'^V • {pSI{p * A/").) 
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In the abscnsc of the last term, equation (|5.32p is a Fokker-Plank equation 

M. = AAi" + /3V.(MV^), (5.34) 

which is known to converge to the stationary solution ha exponentially, where fiA has the mass A := 
||^(-,0)||i and satisfies 

/x™"i = (C-/3^)+ for some OO. (5.35) 



m — 1 



In Therorem 15.131 we will prove for to < 2 — 2/(i, if the initial data is sufficiently less concentrated 
than iJA, then fi{-, r) also converges to the same limit fj,A exponentially as r ^> oo. We begin by defining 
the following mass functions: 



M''(r,r) := A/(r,T;Ai) and X(r,r;/) := / f*h{-,T)dX, 

where M is as given in (j5.ip . h is as given in (|5.33p . and / is an arbitrary function. Note that for 
V = TV, h{-, t) is the delta function for all r, hence M = M. 



Then satisfies the following PDE in the positive set of fi: 

, M(r -r- 1 

.36) 



m ~ 1 dr dr adr ^ (TdT 



We first check that the mass comparison holds for re-scaled equations: 

Proposition 5.11. Let V{x) he given by (A) or (B), and let m < 2 — |. Assume /ii(A, r) is a 
suhsolution and fi2iX,T) is a supersolution of (j5.36p . Further assume that J fii(-,T)dX and J /X2(-,T)(iA 
stay constant for all t>0. Then the mass is ordered for all times, i.e., 

if /ii ( A, 0) -< /i2 ( A, 0) , then we have /ii ( A, r) -< /i2 ( A, r) for all t > 0. 



Proof. Let pi{x,t) be the corresponding re-scaled versions of pLi. Then pi and p2 are respectively a 
subsolution and a supersolution of (|5.2p . The proof then follows from Proposition 15 .31 and from the fact 
that 

M{r, r; p,) = e^^-^^^ Mire^\ e^; p,). 

□ 



Next we state a technical lemma which is used later in the proof of the convergence theorem. The 
proof is in the appendix. 

Lemma 5.12. Let k{t) solve the ODE 

k\t) = Cifc(l - fc") + C2k'^+^e-P\ (5.37) 

where Ci,C2,Q;,/3 are positive constants. Then there exists a constant (5 > such that if < fc(0) < 6, 
then k{t) — > 1 exponentially as t — >■ oo. 
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Now we are ready to prove the main theorem. We wiU first prove it for radiahy decreasing solutions. 

Theorem 5.13. Let V{x) be given by (A) or (B), let I < m < 2 — ^ ^.''^d let /i^ be as given in (|5.35|) . 
Suppose fJ-oW 'is radially decreasing, compactly supported and has mass A. Then there exists a constant 
S > depending on d, m, /ig and V , such that if 

fioW ^ S''^lA{sx), 

then the weak solution fJ.{X, r) of (j5.32p with initial data fiQ exists for all t > 0. Furthermore, M{r, r; /i) 
defined in ()5.ip converges to M{r, r; ^a) exponentially as t oo and uniformly in r. 



Proof. The proof of theorem is analogous to that of Theorem 15.61 we will construct a self-similar 
subsolution (/)(A, r) and supersolution r]{X,T) to (|5.32p . both of which converge to /i^ exponentially. 

Observe that (|5.32p can be written as a transport equation 

V • itiv) = 0, 

where the velocity field v is given by 

V VCm™"^) +(3\ + e(i-")^V(M * (AA * h{y, r)). 

m — 1 

Hence the inward velocity field w(r, t; fi) := —v ■ j|j for the rescaled PDE (|5.32p is 

/ \ m d , „, -I, ^ ,._sM{r,T;^i) 

v{r,T;fi) = 77r^ ')+f3r + e<'' ■ 

m — 1 or (Jdr 

We first construct a subsolution (f>{X,T) with the scaling factor fc(r) to be determined later: 

0(A,r) ■.= k''{T)^lA{k{T)X). 
Since fiA satisfies (|5.35p , the inward velocity field of is then given by 

v{r, r; 0) = (1 - k^^"^-'^+')l3r + e(i-")-M!lIi^. 

Note that the last term of v{r,T;(j)) is always non-negative, and thus v{r,T;(j)) > (1 — fc'''^'"~^^+2)/3r. 
That motivates us to choose fc(r) to be the solution of the following equation 

fc'(T) = - fc'i("-i)+2)^ (5.38) 

with initial data fc(0) sufficiently small such that 0(-,O) -< ^a and </>(•, 0) -< /i(-,0). One can proceed as 
in the proof of Theorem l5.6l to verify is indeed a subsolution. Moreover, it can be easily checked that 
^(t) — > 1 exponentially as r — > oo, hence M(r, r; </>) converges to M(r; //^) exponentially as r — > oo and 
uniformly in r. 

Next we turn to the construction of a supersolution of the form 

7y(A,T) ■.= k'\T)^,A{HT)X). 
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Here the main difBculty comes from the aggregation term, which might cause finite time blow-up of the 
solution. To find an upper bound of the inward velocity field, we first need to control A4{r, r, k'^fiAikX))'- 

JB(0,r) 

< \\AV\\i [ k'^^iA[k\)d\ 

= \\AV\\i f MA(A)dA <C{krf/ad, 

JB{Q.kr) 

where the first inequality is due to Riesz's rearrangement inequality and the fact that fiA is radially 
decreasing. C is some constant that docs not depend on k,r,T. 

The above inequality gives the following upper bound for the inward velocity field of 77: 

v{r, T- 7?) < (1 - fc'^^™-^)+2)^r + Cfc''e(i~")^r. 

Therefore if we let k{t) solve the following ODE 

k'{T) = /3fc(l - fc''('"-i)+2) + C•fc''+le(l-")^ (5.39) 

and choose the initial data fc(0) such that /i(-,0) -< f7(-,0) = fc'^(0)/iA(fc(0)A), then r] would be a 
supersolution to (|5.32p . 

Let us choose fc(0) = 5, where 5 is given in the assumption of this theorem. Due to Lemma [5.12[ 
k{T) — > 1 exponentially as r — >• 00 when 6 is sufficiently small, hence it follows that M(r, r; 77) converges 
to A/(r; ^a) exponentially. 

Since the supersolution rj exists globally, we claim that the weak solution exists globally as well. 
Suppose not: then due to Theorem 4 of jBRBj . ^ has a maximal time interval of existence T*, and 
limT-/iT* ||/^(-, t) lloo = 00. On the other hand. Proposition 15.111 yields that 

r) 77(-, r) for ah r < T*. (5.40) 

Note that Proposition |42] implies that ji is radially decreasing for all t <T*, which gives 

M;r)\\oo < \H;r)\\oo for aU r < T* . (5.41) 

The above inequality implies that IrniT-y^T' ||'y('; t) l|oo — 00, which contradicts the fact that ||'?(-, t) l|oo 
is uniformly bounded for all r. 

Once we have the global existence of fi, Proposition 15.111 vields that 

0(-, r) -< fi{-, t) -< ■/]{■, t) for aU r > 0. 

Since both and rj converge exponentially towards jiA as r — > 00, we can conclude. □ 

The following generalization of Theorem 15.131 will be proved in Section 6. 

Corollary 5.14. Let V{x) he given by (A) or (B), and 1 < m < 2 — ^. For a nonnegative function hq 
in L^{M.'^), define A :~ J jif){X)dX, and let Ha{X) be as given in (|5.35l) . Then the following holds: 
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(a) there exists a small constant 6 > depending on d,m,fiQ and V , such that if 

then the weak solution /.i(A,T) of (|5.32[) with initial data exists for all r > 0. 

(b) Let fiQ is as given in (a) and also is radially symmetric and compactly supported, then M{r,T]^) 
defined in (|5.ip converges to M{r,T; (ia) exponentially as t oo and uniformly in r. 

If we rescalc back to the original space and time variables, Theorem 15.131 immediately yields the 
algebraic convergence of mass function for the solution to 

Corollary 5.15. Let V,m,fi and fiQ be as given in Corollary \5.14\ and let p be given by (j5.3ip . Let 
bl{x,t) be as given in (|5.30p . Then p is a weak solution to (|l.ip . and p vanishes to zero as < — > cxd with 
algebraic decay. In particular if pQ is radially symmetric then 

(a) \M{r,t) — M{r,t;U)\ < Ct^^ , for all r > 0, for some C, 7 depending on po,m,d and V. 

(b) for all p > 1 we have 

where C, 7 depend on p{x,Q),m,d and V . 

Proof. From the proof of Theorem l5.131 we have Wp{ ^^'j^'^ , ^) < e"''"^ for some 7 depending on p{x, 0), 
m, d and F, and the proof is analogous with the proof of CoroUarv 15.81 Now we scale back, and the 
above inequality becomes 



□ 



6 A comparison principle for general solutions and Instant reg- 
ularization in L°° 

In this section we consider general (non-radial) solutions of (jl.ip . Our goal is to prove the following 
result: 

Theorem 6.1. Suppose m > 1. Let V be given by (A) or (B), and let p be the weak solution to (|l.ip 
with initial data p{x,0) ~ Po{x). Let p be the weak solution to (jl.ip with initial data 0) — Pq{x). 
Assume p exists for t £ [0,r), where T may be infinite. Then p*{-,t) -< p{-,t) for all t G [0,r). 

As an application of Theorem 16. 11 we will show that solutions of (|l.ip with initial data in L^ imme- 
diately regularize in (see Proposition 16.61 ) 

The proof of Theorem 16. 1[ which we divide into several subsections follows that of the corresponding 
theorem for solutions of (|1.2p (see Chapter 10 of [V])- The theorem in |V] is proved by taking the 
semi-group approach and applying the Crandall-Liggett Theorem. The challenge lies in the fact that 
our operator in (jl.l[) is not a contraction, in either L^ or L°° . For this reason the proof requires an 
additional approximation of our equation with one with fixed drift: see (|6.4p . 
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6.1 Implicit Time Discretization for PME with drift 



Consider the following equation 

pt = Ap™ + V- (pV$), (6.1) 
where <!> is a function given a priori such that G C(W^ x [0,cxd)), and G C'^iW^) for all t. 

Following the proof in the case of (|1.2p in |V], we approximate (|6.ip via an implicit discrete-time 
scheme. For a small constant h > 0, Ui is recursively defined as the solution of the following equation: 

= AUr + V ■ (C/,V$,), ^ = l,2,... (6.2) 

where ?7o = u{-, 0), <I>i = $(•, i/i). Now define 

p^(.,f) := [/i(.) for (i - l)/i< t < i/i, i = l,2,... (6.3) 

The following result states that our approximation scheme is valid: the proof is in the Appendix. 

Proposition 6.2. Let uq G ^^(IR''; (1 + \x\'^)dx) n L°°(R'^), and let ph be defined by ^^.Then there 
exists a function p G L°°{[Q, oo)\ L^{W'')) such that 

sup ||p(-,i) - Ph(-,t)|lLi(Kd) 
0<t<T 

for any T > 0. Moreover, p coincides with the unique weak solution for (|6.ip . 



6.2 Rearrangement comparison 

For a given function u{x) : M'' — > M, let us define u* as given in (|1.6I) . 

Consider the following equation, where f{x,t) G C([0, oo); L^(M^)) is a given function: 

pt = Ap"' + V- (pV(/*y)). (6.4) 

Theorem 6.3. Suppose m > 1. Let 6e given by (A) or (B), and let p be the weak solution to (j6.4p 
with initial data p(x, 0) = pi:){x). Let p be the weak solution to the symmetrized problem 

Pt = Ap™ + V • (pV(/* * V)), (6.5) 

with initial data p(x, 0) = Pq{x). Then p is radially decreasing, and 

p*{-,t) -< p{-,t) for all t > 0. 

Due to Proposition 16. 2[ to prove Theorem 16.31 it suffices to show the following Proposition; see the 
Appendix for the proof. 

Proposition 6.4. Suppose V is given by (B) and m > 1. Let u £ D (the domain D is defined in 
(|A.lip ) be the weak solution of 

- HAu"" ~ W • (uV(/ *V)) + u = g, (6.6) 

where f,g G L^IW^) are nonnegative. Also, let u £ D be the solution to the symmetrized problem, i.e. 
u solves (j6.6p with /, g replaced by f and g respectively, where f and g are radially decreasing, have the 
same mass as f and g respectively, and satisfy f*~<f and g* < g. Then u* < u. 
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Proof of Theorem [6731 



The radial monotonicity of p can be shown via a similar argument as in Theorem 14.21 in fact the 
argument is easier here since /* * Al/ is a radially decreasing function. 

Next we prove p* < p for all t > 0. Let Ui be the discrete solution for the original problem, and let Vi 
be the discrete solution for the symmetrized problem. Due to Proposition 16.21 it suffices to prove that 
U* -< Vi for alH G N. Here Ui solves 

= A[/r + V ■ ([/.V(/. * V)), (6.7) 
where Uq = u{-,0),fi = f{-,ih), and Vi solves 

Yi^^ = Avr + y-{vy{f*^v)), (6.8) 

where Vo = u*{-,0). Since ~< Vq, by applying Proposition 16.41 inductivelv we can conclude. Lastly 
when V = M , we can use a mollified Newtonian kernel to approximate A/", and the result follows via 
Proposition □ 

Now we are ready to prove our main result: 

Proof for Theorem I6.lt Let us first prove the theorem when V is given by (B), where we have 
global existence of solutions. Let pi{-,t) := p*{-,t) for all t > 0, where p{x,t) is the weak solution 
of (jl.ip with initial data p{x,0) = po{x). For z > 1, we let pi be the weak solution to the following 
equation: 

(p,)t = A(p,)" + V • (p, ■ V(p,_i * V)), (6.9) 
with initial data Pi{x, 0) = p*{x,Q). Observe that pi{-, t) is radially decreasing for all i G N+, i > 0. 

By Theorem 16.31 we have pi -< p^+i for all i G N. Hence we have 

p*{-,t) = pi{-,t)^P2{-,t)^p^{-,t)^..., foralH. (6.10) 

Due to Theorem l3.1l {pi} is locally uniformly continuous in space and time. Hence by the Arzela-Ascoli 
Theorem any subsequence of {pi} locally uniformly converges to a function p along a subsequence. On 
the other hand p is the unique weak solution for (jl.ip with initial data p{x,Q) = Pq{x). This means 
that the whole sequence {pi} locally uniformly converges to p. Now we can conclude due to (|6.10p . 

When V = TV, we can use a mollified Newtonian kernel to approximate A/", and the result follows via 
Proposition 13. 31 □ 

Corollary 6.5. Suppose m > 1. Let V be given by (A) or (B), and let p be the weak solution of p.l[) 
with initial data po{x). Let p be the solution to the symmetrized problem, i.e. p is the weak solution 
to (|1.1|) with initial data Pq{x). Assume p exists for t G [0,r), where T may be infinite. Then for any 
p G (1, oo] we have 

\\p{-,'t)\\Lp{W) < \\p{-,t)\\Lp{W), for all t G [0,r). 



We arc now ready to generalize Theorem 1 5 . 1 31 to non-radial solutions. 

Proof of Corollary I5.14t Let p.{X,T) be the weak solution to (|5.32p with initial data p-oi^). 
Then /i(-,0) meets the assumptions for Theorem 15. 131 which implies the global existence of p. Due to 
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Corollarv l6.5[ \\fi{-, t) |[oo < WfJ-i'i ''") lloo for all t during the existence of /i; hence the uniform boundedness 
of fl yields that fi cannot blow up and thus must exist globally in time. This proves (a). 

Now suppose fiQ is radially symmetric and compactly supported, and fiQ satisfies the assumption 
in (a) such that the corresponding solution fi exists globally in time. In this case we can construct 
subsolution and supersolution as in the proof for Theorem 15.131 to prove (b). □ 

6.3 Instant regularization in L°° 

We finish this section by presenting the following regularization result as a corollary of Theorem 16. II It 
says that for initial data po S L^{M.'^; (1 + |xp)(ia;) fl i°°(M'^), no matter how large the L°° norm of po 
is, ||p(-,i)||oo will always be bounded by t~" for some short time. 

Proposition 6.6. Let V be given by (A) or (B), and let m > 2 — 2/d. Let p{x,t) be the weak solution for 
([LT]) . with initial data po G L^{R'^; {l + \x\'^)dx)nL°°{M.'^). Let us denote A ^ \\po\\i and a d(m-i)+2 - 
Then there exists c ~ c{m, d, A, V) and to = (2c)"'^/" > such that we have t) G L°°(M'') with 

l|p(-,i)||Loo(Rci) < c{m,d,A,V)t-°' for all < t < to- 

Proof. By Corollarv l6.51 it suffices to prove the inequality when po is radially symmetric. Also, in this 
proof we denote c(m, d, A, V) to be all constants which only depend on m, d, A, V. 

Let PA be the radial stationary solution of (|l.ip with mass A. Note that pA is radially decreasing, 
and thus pa(0) > 0. Since uq is a radial function in we can scale pA to make it more concentrated 
than Uq, i.e. we choose < a < 1 to be sufficiently small such that 

Uq -< a^'^pAia^^x). 
As in the proof of Theorem 15.61 let us define 

7?(r,i) -.^k'^it) pA{Ht)r), 
where k{t) solves the following ODE with initial data fc(0) = a^^: 

k'{t) ^ c{m,d,A,V)k'^+\l - fc'i(™-2+2/d))_ 
Here c(m, d. A, V) corresponds to C2 in the proof for Theorem 15.61 It was shown in the proof that 

p{-,t) <T]{-,t) forant>0, 

which in particular yields that 

IkK-.i)llL-(K'i) > l|p(-,OllL~(Rrf) for all t > 0. 

Observe that, by definition, 

Ht) := H;t)\\L^m = k''{t)pAiO) = c{m,d,A,V)k''it). 

Therefore to prove our proposition it is enough to show 

h{t) < f{t) := c{m, d, A, for aU /i(0) > and t G [0, to] , 



28 



where to is chosen such that f{t) > 2. Note that h{t) solves 

h'{t) = c{m,d,A,V)k'^-'^k' 

^ c{m,d,A,V)h\l~h'"~^+^/''). 

In particular when h(t) > 2, h satisfies the following inequality 

h'{t) < -c{m,d,A,V)K"'+'^''^. 

Since f{t) solves the above ODE with equality, we obtain h{t) < /(t) for < t < tg. Now we are 
done. □ 



A Appendix 

A.l Proof of existence for p as given in Proposition 2.1 

Here we will show the existence of a global niinimizer for the free energy functional given in (|1.5p . First 
note that the kernels V given in (A) and (B) belong to M'^^ , where denotes the weak space. 

Our proof is based on a theorem of Lions in [L] : 
Theorem A.l ([L]). Suppose f G MP{M.'^), f > and consider the problem 



Ix = inf ( / u™d.T - - 

u£Kx L /iRd m — 1 2 



where 

Kx = \u e i«(R'*) n L^{W^), M > a.e., / udx = a| with q = 'EJlI^ 
Then there exists a minimizer of problem (Ix) if and only if the following holds: 

Ix<Ic.+h-c., Vae (0,A). (A.l) 

Proposition A. 2 ([L|). Suppose there exists some a G (0,d) such that f{tx) > t^"f{x) for all t>l. 
Then \A.1\) holds if and only if 

Ix < 0, for all A > 0. (A.2) 

For the rest of this subsection, we will verify that Proposition IA.2I applies to our kernels. 

Proposition A.3. Let f ^ V given by either ([0|) or pr4| . Then f{tx) > t-°'f{x) with t > 1 and 
a^d-2. 

Proof. When V is given by (A) the proof is straightforward, so suppose V is given by (B). Then 
/ = —A/" * h. Since h can be approximated by a sum of indicator functions, it suffices to prove the 
proposition for / = —A/"* XB{o.r)j where x is the indicator function. In this case we have 

1 „ 1 



-r 



2 —\x\^ for |a;| < r, 



/(.)=^2(d-2) 2d' (^3) 
___,,d|2.|-d+2 foj. |x| > r, 

which finishes the proof. □ 
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Proposition A. 4. Let f be as in Lemma ] A. 3\ Suppose to > 2 — | and let us define 

u = —r^ 

where Cd is the volume of the unit ball in M.'^ and R is a constant to be chosen later. If R is sufficiently 
large, we have 

E{u):= ( u^dx-i / u{u*f)dx<Q, 



TO - 1 2 



and thus I\ < E{u) < 0. 



Proof. First note that we have 

V cR'^ J 



1 ..m^^^ f 1 f^XBiO,R)\"'^ 



TO — 1 /md m — 1 V cR"^ 



On the other hand, /g^p _r/2)(~^)'^^ ~ i?^ if i? is sufficiently large: this implies 

/jgd u{u * {-V))dx > A^i?^'*+^. Since m > 2 - | we have 



Elu) = / u'^dx - u(u* (~V))dx < 0. 

' m — 1 ' 



□ 



A. 2 Proof of Lemma 14.11 

Observe that AV is nonnegative and radially decreasing, and thus it can be approximated in L^{M.'^) D 
L°°{M.'^) by the sum of bump functions of the form cxB{o,r), where c > 0. By linearity of convolution, 
it suffices to prove that for each bump function X-B(o,r)j where r is any positive real number, we have 

{U * XB{0,r)){bl) - (m * Xi3(0,r))(ai) < || XS(0,r) II 1 ) " "(oi)). (A.4) 

Observe that 



{u * XB{o,r)){bi) - (u * XB{o,r)){ai) = / u{x)dx-l u{x)dx (A.5) 

JB{bi,r) JB{ai.r) 



u{x)dx — I u{x)dx. (A. 6) 

JflA 

Here ftA := B{ai,r)\B{bi,r) and ^Is := B{bi,r)\B{ai,r) (see Figure gj). 

Note that fi^ and fis are symmetric about the hyperplane H — {x : xi ~ ^^^}- For any x G fi^, 
use f{x) to denote the reflection point of x with respect to H. Then we have 



u{x)dx — / u{x)dx = I {u{f{x)) — u{x))dx. 
Since |a;| < |/(a;)| for x G Q,a, we can use the assumption ()4.ip to obtain 

{u{f{x))-u{x))dx < I {u{b)-u{a))dx<\B{Q,r)\{u{b)-u{a)), 
which completes the proof. 
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_a+/3 
Xi — 2 



Figure 4: The domains ft a and fi^ 
A. 3 Proof of Lemma [5.121 

Proof. When < fc < 1, the right hand side of (|5.37p is bounded above by (Ci + C2)fc. Hence if the 
initial data satisfies < fc(0) < 1, the inequahty k{t) < fc(0)e''-^i+'^^-'* wih hold until k reaches 1. In 
other words, k{t) is guaranteed to be smaller than 1 until time ti :~ — p'^^, ■ 

Now if we choose fc(0) to be sufficiently small such that < k{Q) < S, where 

C1+C2 

then ti would be sufficiently large such that 

- 2 

We claim g{t) := 1 + e^'^^*~^^^ is a supersolution of (|5.37p for t > ti, where e := min{/3, ^Cia}. It is 
obvious that g{ti) > 1 > k{ti), so it suffices to show 

5'W >C^iff(l-5") + C2/+'e-'5* fori>ti. (A.7) 

By definition of g, we have 

RHSofdMl < -Ciae-'(*-*i) +C22'*+ie-''*ie-''(*-*i) (A.8) 

< -iCiae-'(*-*i) (A.9) 

< -ee-'(*-*i) =:LHSof(IMl- (A.IO) 

Therefore k{t) < 1 + e-<=(*-*i) for ah t > h. 

To obtain the corresponding lower bound for k{t), note that the last term of (|5.37p is non- negative. 
Therefore if g solves g'{t) = Cig{l — g'^) and g(0) = fc(0), then \g{t) — 1| ^ e"'-^^"*. Comparison between 
these two ODEs yields k{t) > g(t) for all t > 0, which implies k{t) > 1 — Ce"*-^^"*. Now we can conclude 
that there exists C depending on Ci,C2,a,/3 and fc(0) such that 

\k{t) - 1| < Ce"'*. 

□ 
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A. 4 Proof of Proposition 16.21 



The proof of ProDOsition l6.2l is an application of the Crandall-Liggett Theorem f |CLj . also see Theorem 
10.16 in [V]). Let us consider the following domain: 



Wi^;^(]R''), Am™ g L^W^), |Vii™| G M'^/^'^-^\R'^)y (A.ll) 



Here the Marcinkicwicz space Mp{R'^), 1 < p < cxo. is defined as the set of / G Lj^^CM.'^) such that 



\fix)\dx < C\K\^P-^^/P, 



IK 

for all subsets K of finite measure. The minimal C in the above inequality gives a norm in this space, 



) = sup{meas(i^)-(P-i)/'' J \f\dx : K C M'',meas(ii") > o|. 



\\J llMP(K<i) — 

A parallel argument as in Theorem 2.1 of |BECj yields the existence of solutions for the discrctized 
equation. 

Lemma A. 5 (Existence). Let d > 3 and let uq G L"'^(R'^),$ G C^(R'*). Then there exists a unique 
weak solution u ^ D of the following equation: 

^^-^ Am'" + V-(uV4>). (A.12) 

The proof of the next lemma is parallel to that of Prop 3.5 in [V] for (|1.2p . 

Lemma A. 6 (L^ contraction). Let <I> G (R'^), Uq, G L^iR"^) and let Ui,U2 E D be the weak solutions 
to the degenerate elliptic equation 

^^^^^ = A(u,)'" + V-(u,V$), i = l,2. (A.13) 
h 

Then ui and U2 satisfy 

\\ui - U2||Li(Rd) < ||uoi - M02||li(B<')- (A. 14) 

Proof of Proposition [572] 

Proof. Let D be defined above, and define the nonlinear operator A: D ^ L^(M'*) by 

A{u) = -Au™ - V(uV$), 
Then Lcmma lA.51 and Lcmma lA. 61 yield that for any /i > 0, there is a unique solution uin D solving 

hA{u) + u = f. 

Moreover the map / m is a contraction in L^(R''). Now arguing as in |V], the Crandall-Liggett 
Theorem yields the conclusion. □ 
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A. 5 Proof of Proposition 16.41 



The proof of Proposition |6]4] is parallel to that of Theorem 11.7 in [V] for (|1.2p . First we state a lemma 
which deals with the extra convolution term. 

Lemma A. 7. Let V be given by (B). Let f G L^M.'^) and e W^'^^iR'^) be non-negative functions. 
Then for any non-negative number a, b, we have 



V(/ * (-V^)) ■ V0 < / if* *AV){max{(t)*,b}~a), (A.15) 

{a<4><b} J{4>*>a} 

where the equality is achieved if f, (j) are both radially decreasing. 
Proof. Let us define ry : M'^ — M by 

f b if > b, 

r]{x) := < (j){x) — a if a < (j){x) < b, 
[q if0(.T)<a. 

Then T]{x) £ W^'°°{R'^), = Vt? in {a < <l){x) < b}, and Vt? = in R'^\{a < (j){x) < b}. Therefore 
LHS of (IXTSI) = / V(/ * {-V)) ■ V77 



< {f*^AV)rf= (/* *Al^)(max{0*,6}-a), 

jRd J{0'>a} 

where the inequality comes from Riesz's rearrangement inequality. Note that we obtain an equality if 
/ = /* and r] = T]*. Hence the lemma is proved. □ 

The following lemma corresponds to Theorem 17.5 in [V] . 

Lemma A. 8. Let V be given by (B). Let f,f and g be non-negative radially decreasing functions in 
L^{R'^), where f < f ■ Let h > 0, and let vi,V2 £ D be two non-negative radial decreasing functions. 
Assume vi and V2 satisfies 

- hA{v,r - hV ■ (i-i V(/ * V)) +vi-<g, (A.16) 

- hA{v2r - hV ■ iv2y{f * V)) -\- V2 = g. (A.17) 

Then we have vi -< V2. 

Proof. Let Ui := v™ and define m ;= ui — M2, w := wi — W2, Mr) :— /g^p v{x)dx. Our goal is to show 
A{r) < for all r > 0. 

Subtracting (jA.16|) from (jA.17|) . and integrating the quantity in B{0,r) yields that 

/ -hAudx - h(vi{r) f*AVdx-V2{r) f * AVdx) + A{r) < 0, (A.18) 

J B(0,r) ^ J B(0,r) J B{0,r) ' 
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which can be written as 

-hcdr^-^u'{r) -hv{r) ( f * AVdx - hv2{r) [ {f - f) * AVdx + A{r) < 0. (A.19) 

J B(0,r) J B{0,r) 

(Here u'ir) exists due to the fact that G I? for i = 1, 2, which imphes that Au is in L^.) Since we 
assume / ^ /, it follows that Jg^j, r)^{f ~ f) * AV)dx < for all r > 0. Therefore 

- hcdr'^-^u'{r) - hv{r) / * AF + A(r) < for all r > 0. (A.20) 

JB(0,r) 

Note that since Ui and Vi both vanish at infinity, from (|A.20p it follows that limr_>oo A{r) < 0. Hence if 
A{r) is positive somewhere, it achieves its positive maximum at some point ro > 0. At r = vq we have 
(ro) = A'(ro) = 0, and (|A.20p becomes 



V 



which means U2 — ui is strictly increasing at vq: hence V2 — vi will also be strictly positive in (ro, ro + e) 
for some small e, which implies A{rQ + e) > A(ro). This contradicts our assumption that A(r) achieves 
its maximum at ro. Therefore A{r) < for all r, which means U2 -< ui- □ 



Proof of Proposition [674t The proof is parallel to that of Theorem 11.7 as in |V]. For any test 
function (jj e Wo'°°{R'^), we have 

h[ Vu"-V0 + /i/ u\/[f *V)-y(f)+ I u(f)^ I gcj), (A.21) 

where (j) G Wq'°°{W'') is any test function. Now let us take (j>{x) :— (w™(x) — i)+ where t > 0, and 
differentiate the equation with respect to t. Then wc have: 

^h{^ f |vu"f / ^v(/*y).v(ri'"+^))+ / [ g. (A.22) 

h h h li 

Following the proof of Theorem 17.7 in [V], one can check that 

h < I /iA((u*)™) (with equality if u EE li*). 



^3 = / U*, 



/4 < sup 9* ^ 9*- 

\n\=\{u^>t}\Jn J{{u')™>t} 

It remains to examine I2 ■ Using Lemma IA.71 it follows that 

I2 ^ h hm - I ^^V(/ * {-V)) ■ V(u"+i) 



< /iliminfi / ^^(r *Al/)(max{w'"+\(t + e)i+^}-ti+^)+ 

e J{t<(«-)'"<t+e} '71+ 1 



/li ™ / /* * Ay. 
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Plugging in the four inequalities into (|A.22p . the following inequality holds for alH > 0: 



hA{{u*y') ~ ht- / f**AV+ u* < g*. (A.23) 

Since f > is arbitrary, the above inequality implies 

- /iA((u*)™) - W ■ (w*V(/* * Vj) +u* < g*. (A.24) 

On the other hand, by assumption, u solves 

- /iA(u™) - hV ■ {uV{f *V)) + u^g, (A.25) 

where f >~ f* and g > g* . Note that u <E D implies u* G D. So we can apply Lemma [A. 81 and get 

u* ~<u. □ 
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